Elliptic flow in ]9]9-collisions at the LHC 



S.M. Troshin, N.E. Tyurin 

Institute for High Energy Physics, 
Protvino, Moscow Region, 142281, Russia 



Abstract 

We consider collective effects in pp-coUisions at the LHC energies re- 
lated to presence of the large orbital angular momentum in the initial state 
and role of this orbital momentum in the elliptic flow behavior. 
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Introduction 



There were several collective effects in nuclear collisions observed at RHIC and 
various experimental observables relevant to their studies have been measured. In 
particular, there was observed ridge structure in the near- side two-particle corre- 
lation function in the secondary particles production. Similar structure in the two- 
particle correlation function was observed by the CMS Collaboration [HI. This is 
not an expected result. Commonly, pj9-collisions were treated as a kind of "ele- 
mentary" ones or the reference process for detecting deconfined phase formation 
in AA-coUisions. It becomes evident that such a view should be corrected. 

The experimental observations of RHIC and LHC on the ridge in the two- 
particle correlation functions demonstrate that the hadronic matter is strongly cor- 
related and reveals high degrees of coherence when it is well beyond the critical 
values of density and temperature. The other anisotropic flow measurements at 
RHIC also demonstrated such a collective behavior. 

In this note we discuss the elliptic flow in proton collisions on the base of 
approach [|2l which was applied for discussion of the directed flow and the ridge 
in proton collisions f3]. This particular approach is based on the non-perturbative 
hadron dynamics. It should be noted that the values of transverse momenta in the 
relevant experimental data are not too large. Despite that there are considerations 
of the experimental results rendering to the perturbative QCD and parton radiation 
in the medium. 



1 Elliptic flow in peripheral proton collisions 

We consider in this section peripheral hadronic collisions . These collisions sup- 
pose that the impact parameter is different from zero. The mechanism of selecting 
such events in proton collisions will be considered further for the LHC energies. 
Here we adopt that impact parameter is non-zero and the reaction plane can there- 
fore be determined by the cumulant method. 

There are several experimental probes of collective dynamics. A most widely 
discussed one is the elliptic flow 

V2ip±) = (cos(20))p, = {^^^^), (1) 

p± 

which is the second Fourier moment of the azimuthal momentum distribution of 
the particles with a fixed value of pj_. The azimuthal angle is an angle of the 
detected particle with respect to the reaction plane, which is spanned by the col- 
lision axis z and the impact parameter vector b. The impact parameter vector b 
is directed along the x axis. Averaging is taken over large number of the events. 
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Elliptic flow can be expressed in covariant form in terms of the impact parameter 
and transverse momentum correlations as follows 



^2 P± = 2 - 2 ' (2) 

PI PI 

where b = h/h. To get some hints on the possible behavior of the elliptic flow 
in proton collisions, it is useful to recollect what is known on this observable 
from nuclear collisions experiments. Integrated elliptic flow V2 at high energies is 
positive and increases with ^/sNN. The differential elliptic flow f2(p±) increases 
with p± at small values of transverse momenta, then it becomes flatten in the 
region of the intermediate transverse momenta and decreases at large pj_. 

It is also useful to apply a simple geometrical ideas which imply existence 
of the elliptic flow in hadronic reactions. Geometrical notions for description of 
multiparticle production in hadronic reactions were used by many authors, e.g by 
Chou and Yang in lIH. In the peripheral hadronic collisions the overlap region 
has different sizes along the x and y directions. According to the uncertainty 
principle we can estimate the value of p,j. as 1/Ax and correspondingly Py ~ 
1 / Ay where Ax and Ay characterize the size of the region where the particle 
originate from. Taking Ax ~ Rx and Ay ~ Ry, where Rx and Ry characterize 
the sizes of the almond-like overlap region in transverse plane, we can easily 
obtain proportionality of V2 in collisions with fixed initial impact parameter to the 
eccentricity of the overlap region, i.e. 

f?2 d2 

^2 ~ ^ ~ (3) 

R^ -\- Ry 

The presence of correlations of impact parameter vector b and px in hadron 
interactions follows also from the relation between impact parameters in the mul- 
tiparticle production: 

h = ^Xihi. (4) 

i 

Here Xj stand for Feynman xp of i-th particle, the impact parameters bj are con- 
jugated to the transverse momenta Pi j_. 

The above considerations are based on the uncertainty principle and angular 
momentum conservation, but they do not preclude the existence of the dynamical 
description, which will be discussed in the next section 
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2 Peripheral inelastic proton-proton interactions at 
the LHC energies 

The particle production mechanism proposed in the model [12] [Si takes into ac- 
count the geometry of the overlap region and dynamical properties of the transient 
state in hadron interaction. This picture assumes deconfinement at the initial stage 
of interaction. The transient state here appears as a rotating medium of massive 
quarks and pions which hadronize and form multiparticle final state at the freeze- 
out stage. Essential point for this rotation is the presence of a non-zero impact 
parameter in the collision. 

In this section we consider unitarity saturation as a dynamical mechanism 
leading to the peripheral nature of inelastic collisions at the LHC energies. We 
also discuss the limiting energy dependence of the anisotropic flows resulting from 
unitarity saturation in central coUisons (6 = 0). 

In the untarization approach ([/-matrix) the elastic scattering matrix in the 
impact parameter representation has the form: 

l + tUis,b) 

S{s,b) = - , (5) 

1 — iU[s, b) 

where S{s, b) = 1 + 2z/(s, b) and U{s, b) is the generalized reaction matrix, which 
is considered to be an input dynamical quantity similar to the eikonal function. 
Unitarity equation rewritten at high energies for the elastic amplitude /(s, b) has 
the form 

Im/(s, b) = hei{s, b) + hinei{s, b) (6) 
where the inelastic overlap function 

is the sum of all inelastic channel contributions. Inelastic overlap function is re- 
lated to U{s, b) according to Eqs. (jSj and Q as follows 

, ImUis.b) 

^-^'^^'^) = ii -.[/(., 6)r 

this function is the probability distribution over impact parameter of the inelastic 
interactions: 

Im[/(s, b) 
It should be noted that 

\mU{s,b)=Y,Un{s,b). (9) 

n>3 



(yinel{s) =871 I bdb- T777^-TVm- (8) 
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where ?7„(s, b) is a Fourier-Bessel transform of the function f7„(s, t), the form of 
this function can be found in [|5]|. 

Since the form of the impact parameter dependence of the function function 
hinei{s, h) evolves from a central to the peripheral one at the LHC energies [6J and 
this function tends to zero at 6 = and s — )■ oo, the mean multiplicity 

/ \/ ^ bdb{n){s,b)hineiis,b) 
jg bdbhinei{s,b) 

obtains a main input from the collisions with non-zero impact parameters. One 
can suggest therefore that the events with average and higher multiplicity at the 
LHC energy ^/s = 7 TeV correspond to the peripheral hadron collisions [6 1. Thus, 
at the LHC energy ^/s = 7 TeV there is a dynamical selection of the peripheral 
region in impact parameter space responsible for the inelastic processes. In the 
nuclear reactions similar selection of the peripheral collisions is provided by the 
relevant experimental adjustments. 

Of course, the standard inclusive cross-section for unpolarized particles being 
integrated over impact parameter b, does not depend on the azimuthal angle of 
the detected particle transverse momentum. It can be written with account for 
s-channel unitarity in the following form 

da /'°°, „ I(s,b,u) 

;r = 8^ / bdb \\ i (10) 

du Jq \l-tU{s,b}\^ 

When the impact parameter vector b and transverse momentum of the 
detected particle are fixed the function / does depend on the azimuthal angle (p 
between vectors b and p^. It should be noted that the impact parameter b is the 
variable conjugated to the transferred momentum q = pjj — Pa between two inci- 
dent channels which describe production processes of the same final multiparticle 
state. The dependence on the azimuthal angle can be written in explicit form 
through the Fourier series expansion 

^ oo 

^(■s,b,?/,p_L) = —Io{s,b,y,p±)[l + y^2vnis,b,y,pj_) cosn0]. (11) 

n=l 

The function /o(s, b,^) satisfies to the following sum rule 

j Io{s,b,y,p±)p±dp±dy = {n){s,b)lmU{s,b), (12) 

where n(s,b) is the mean multiplicity depending on impact parameter. Thus, the 
bare anisotropic flow w„(s, b, y, pj_) is related to the measured flow Vn as follows 

Vn{s,b,y,pA_) = w{s,b)vn{s,b,y,pA_). 
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where the function w{s,b) is 



w{s,b) = \l-iU{s,b)\-^. 

In the above formulas the variable y denotes rapidity, i.e. y = sinh~^(p/m), 
where p is a longitudinal momentum. Thus, we can see that unitarity corrections 
are mostly important at small impact parameters, i.e. they modify anisotropic 
flows at small centralities, while peripheral collisions are almost not affected by 
unitarity. The following limiting behavior of f „ at 6 = can easily be obtained: 

Vn{s,b = 0,y,p±) 

at s — )■ oo since U{s,b = 0) — oo in this limit. 

In conclusion, one should note that inelastic events at the LHC energies come 
mainly from the non-zero values of the impact parameter. It is the result of ap- 
proaching to unitarity saturation. 

3 Coherent rotation of transient matter and elliptic 
flow 

Before turning on discussion of anisotropic flows, it is useful to note that the 
energy dependence of the average transverse momentum was considered in [|5l. 
The mechanism is the following: the nonzero orbital angular momentum leads to 
coherent rotation of a transient matter located in the overlap region as a whole 
in the x2-plane . This rotation is similar to the rotation of a liquid where strong 
correlations between particles momenta exist. 

In this mechanism of hadron production, the valence constituent quarks ex- 
cite a part of the cloud of virtual massive quarks and those quarks subsequently 
hadronize and form the multiparticle final state. 

Here we use this mechanism for evaluation of the elliptic flow in pp-interactions. 
It was suggested in [5J that the rotation of transient matter will affect the average 
transverse momentum of secondary hadrons produced in proton-proton collisions. 
A coherent rotation of transient liquid-like state results in a contribution to the 
transverse momentum in the form: 

ApT = nLis,b), (13) 

where L{s,b) is the orbital angular momentum and k is a constant which has a 
dimension of inverse length. Note that the power-like dependence of the average 
transverse momentum at high energies has been obtained in good agreement with 
experimental data [5j . 
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Going further, one should note that, in fact, the rotation gives a contribution 
to the x-component of the transverse momentum and does not contribute to the 
^/-component of the transverse momentum, i.e. 



Ap^ = kL{s, b) 



(14) 



while 



Apy = 0. 



(15) 



Thus, assuming that px = Po + Apx (Ap^ ^ po) and py = po, the effect of the 
rotating transient matter contribution into the elliptic flow can be calculated by 
analogy with the calculation of the average transverse momentum performed in 
[|5||. The resulting integrated elliptic flow increases with energy 



where 5c = 0.207. It should be noted here that transient matter consists of virtual 
constituent quarks strongly interacting by pion (Goldstone bosons) exchanges . 

The incoming constituent quark has a finite geometrical size determined by 
the radius rg and interaction radius Rq (Rq > tq). Former one is determined by 
the chiral symmetry breaking mechanism and the latter one — by the confinement 
radius. Meanwhile, it is natural to suppose on the base of the uncertainty principle 
that size of the region where the virtual massive quark Q is knocked out from 
the cloud is determined by its transverse momentum, i.e. R ~ However, 
it is evident that R cannot be larger than the interaction radius of the valence 
constituent quark Rq. It is also clear that R should not be less than the geometrical 
size of the valence constituent quark tq for this mechanism be a working one. 
When R becomes less than rq, this constituent quark mechanism does not work 
anymore and one should expect vanishing collective effects in the relevant region 
of the transverse momentum. 

The value of the quark interaction radius was obtained under analysis of the 
elastic scattering [|71 and it has the following dependence on its mass 



where ^ — 2 and therefore Rq ~ 1 /m, while the geometrical radius of quark 
Tq is about 0.2 fm. It should be noted the region, which is responsible for the 
small-p^ hadron production, has large transverse dimension and the incoming 
constituent quark excites the rotating cloud of quarks with different values and di- 
rections of their momenta in that case. Effect of rotation will therefore be smeared 
off over the volume V/j and then one should expect that {Apx)vf, — 0. Thus, 



f 2 OC S 



(16) 



v^{p±) = {V2)v, 



R 



~ 



(17) 
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at small p^. When we proceed to the region of higher values of px, the radius 
R is decreasing and the effect of rotation becomes more and more prominent, 
incoming valence quark excites now the region where most of the quarks move 
coherently, in the same direction, with approximately the same velocity. The mean 
value {Apx)vji ^ ^ ^^'^ 

v^{p±) = {v^Wj, > (18) 

and it increases with p±. The increase of f ^ with will disappear at _R = rq, i.e. 
atp± > l/rq, and saturation will take place. The value of transverse momentum 
where the flattening starts is about 1 GeV/c for rq ~ 0.2 fm. At very large 
transverse momenta the constituent quark picture would not be valid and elliptic 
flow vanishes as it was already mentioned. 

We discussed elliptic flow for the constituent quarks. Predictions for the el- 
liptic flow for the particular hadrons depends on the supposed mechanism of 
hadronization. For the region of the intermediate values of p± the constituent 
quark coalescence mechanism flU |9l would be dominating one. In that case the 
values for the hadron elliptic flow can be obtained from the constituent quark one 
by the replacement V2 — t- nyV2 and p± — )■ p^/ny, where ny is the number of 
constituent quarks in the produced hadron. 

Typical qualitative dependence of elliptic flow in pp-coUisions in this approach 
is presented in Fig. 1 




P_T 



Figure 1 : Qualitative dependence of the elliptic flow f 2 on transverse momentum 
in pp-coUisions. 

The centrality dependence of the elliptic flow is determined by dependence of 
the orbital angular momentum L on the impact parameter, i.e. it should be de- 
creasing towards high and low centralities. Decrease toward high centralities is 
evident since no overlap of hadrons should occur at high enough impact param- 
eters. Decrease of V2 toward lower centralities is specific prediction of the pro- 
posed mechanism based on rotation since central collisions with smaller impact 
parameters would lead to slower rotation or its complete absence in the head-on 
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collisions. Qualitative dependence of the elliptic flow on the impact parameter 
is similar to the curve depicted in Fig.l where variable b being used instead of 
transverse momentum. 

Conclusion 

We considered the elliptic flow in the proton-proton interactions at the LHC ener- 
gies in the particular nonperturbative approach, where the origin of elliptic flow is 
associated with the effects of rotation of quark-pion transient matter. This mech- 
anism of anisotropic flow might be a leading one in hadron coUisons, since those 
have smaller geometrical extension and the probabilty of hydrodynamical gener- 
ation of elliptic flow is lower compared to the collisions of nuclei. 

At such high energies (LHC) the presence of reflective scattering mode serves 
as a trigger for the mainly peripheral nature of the inelastic interactions and the 
presence of large orbital angular momenta. Since correlations are maximal in the 
rotation plane, a narrow ridge should be observed in the two-particle correlation 
function [3J. 

In conclusion, we would like to remark that performing studies of the multi- 
particle productions in pp-coUisions at the LHC under the scope of the searches 
for the possible existence of the rotation effects would be of a significant interest. 
Such effects would be absent when the genuine quark-gluon plasma (gas of free 
quarks and gluons) being formed at the LHC energies because of all collective 
effects and the anisotropic flows, in particular, should dissapear in this case. 
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